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Ttiin  report  relates  to  utility  of  imaging  radars  for  arctic  surface  guidance.  The 
program  involves  field  measurements  of  backscatter  from  snow,  ice  and  sea  ice  for  a 
wide  variety  of  meteorological  and  surface  conditions.  Laboratory  measurements  of 
the  complex  dielectric  constant  are  presented  also,  to  enable  the  separation  of  the 
roughness  contribution  to  the  backscattei;.  During  the  first  seven  months  of  the  pro¬ 
gram,  the  dielectric  properties  of  sea  ice  were  measured  in  the  laboratory  from  26  tc 
40  GHz,  over  a  representative  range  of  temperatures,  salinities,  and  orientation  As 
reported  in  the  detailed  Appendix,  which  is  a  Master's  thesis, >the  real  part  of  the 
dielectric  constant  is  close  to  3  for  most  conditions,  but  the  Tbss  tangent  can  vary 
over  three  orders  of  magnitude,  depending  upon  the  temperature  and  salinity.  Under 
most  conditions,  an  incident  wave  of  1  cm  wavelength  would  penetrate  less  than  one 
meter  into  the  sea  ice.  Strong  distinction  of  freshly-formed  sea  ice  from  old  sea 
ice  appears  likely.  Backscatter  from  sea  ice  is  currently  being  measured.  Distinc¬ 
tion  of  rough  and  smooth  surfaces  is  possible,  based  upon  preliminary  interpretation 
of  the  data.  continuation  of  laboratory  and  field  measurements  is  recommended,  to 
include  arctic  conditions  throughout  the  year. 
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SUMMARY 


Defense  operations  in  the  Arctic  must  be  capable  of  action  during 
the  conditions  of  poor  visibility  usually  encountered.  Imaging  radars 
can  provide  such  guidance,  but  the  returns  must  be  easily  related  to 
the  topography  of  the  snow-covered  land  and  sea  ice.  In  such  interpre¬ 
tation,  the  relative  contributions  of  the  dielectric  properties  and 
of  the  surface  roughness  should  be  separated.  For  resolution  cell  sizes 
of  the  order  of  a  few  meters,  and  reasonable  antenna  dimensions,  wavelengths 
in  the  1  cm  range  or  smaller  are  necessary.  The  present  program  involves 
field  measurements  of  backscatter  from  snow,  ice,  and  sea  ice  for  a  wide 
variety  of  meteorological  and  surface  conditions.  Laboratory  measure¬ 
ments  of  the  complex  dielectric  constant  are  also  included,  to  enable 
the  separation  of  the  roughness  contribution  to  the  backscatter. 

During  the  first  seven  months  of  the  program,  the  dielectric 
properties  of  sea  ice  were  measured  in  the  laboratory  from  26  -  40  GHz, 
over  a  representative  range  of  temperatures,  salinities,  and  orientation. 

As  reported  in  the  detailed  Appendix,  the  real  part  of  the  dielectric 
constant  is  close  to  3  for  most  conditions,  but  the  loss  tangent  can 
vary  over  three  orders  of  magnitude,  depending  upon  the  temperature 
and  salinity.  Under  most  conditions,  an  incident  wave  of  1  cm  wave¬ 
length  would  penetrate  less  than  one  meter  into  the  sea  ice.  Strong 
distinction  of  freshly-formed  sea  ice  from  old  sea  ice  appears  likely. 
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During  the  reporting  period,  a  fixture  for  field  testing  was  de¬ 
signed,  fabricated,  and  tested.  It  made  possible  the  backscatter 
measurements  on  snow-covered  land  which  were  obtained  and  which  are 
being  analyzed.  Backscatter  from  sea  ice  is  currently  being  measured. 
Distinction  of  rough  and  smooth  surfaces  is  certainly  possible,  based 
upon  preliminary  interpretation  of  the  data.  A  continuation  of  labora¬ 
tory  and  field  measurements  is  recommended,  to  include  Arctic  conditions 
throughout  the  year.  Further  comparisons  of  laboratory  and  field  data 
are  also  recommended. 
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I.  INTRODUCTION 


For  many  decades,  man's  operations  in  the  Arctic  have  been 
hindered  by  the  extreme  cold  and  lack  of  visibility.  Immersed 
in  darkness  for  several  months  of  the  year,  the  Arctic  ice  cap 
presents  deadly  obstacles  to  surface  operations,  in  the  form  of 
irregular  pressure  ridges  of  ice  and  open-water  leads.  Fog  and 
blowing  snow  often  obscure  this  topography.  Defense  operations, 
which  must  be  capable  of  action  under  all  weather  conditions,  can 
still  retain  surface  mobility  if  equipped  with  a  high-resolution 
imaging  radar  system  which  can  penetrate  through  the  fog  and 
blowing  snow,  giving  a  radar  map  of  the  snow  and  ice  features. 

In  the  interpretation  of  the  image  information  contained  in  the 
radar  return  signal,  one  would  like  to  separate  out  the  effects  of 
surface  roughness  and  geometrical  irregularities  from  changes  in 
the  dielectric  properties  of  the  reflecting  surface.  Preliminary 
consideration  of  required  distances,  resolution  cell  sizes,  and 
limits  on  the  physical  size  of  antennas  lead  quickly  to  the  con¬ 
clusion  that  millimeter  wavelengths  are  necessary.  Very  little 
is  known  about  the  millimeter  wave  reflectivity  of  snow,  ice,  and 
sea  ice.  The  present  engineering  research  program  is  directed 
towards  the  measurement  of  the  backscatter  from  actual  Arctic 
terrain;  the  laboratory  measurement  of  the  dielectric  constant  and 
loss  tangent  of  snow,  ice,  and  sea  ice;  and  the  correlation  of  the 
field  and  laboratory  measurements.  The  surface  roughness  will  be 
separated  out  as  an  independent  effect,  and  finally,  suggestions 
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will  be  made  for  the  design  of  obstacle  avoidance  systems  which  can 
perform  this  separation  on  a  continuous,  on-line  basis. 

During  the  first  seven  months  of  the  program,  laboratory  measure¬ 
ments  have  been  made  of  the  dielectric  properties  of  sea  ice  in  the 
frequency  range  of  26  -  40  GHz.  Preparations  for  field  measurements  during 
the  winter  were  completed,  and  field  measurements  of  snow-covered  land 
were  begun.  In  the  Appendix,  a  detailed  report  of  the  laboratory 
measurements  on  sea  ice  is  presented.  The  following  sections  describe 
some  of  the  major  findings  and  their  implications. 

II.  DIELECTRIC  PROPERTIES  AND  RADAR 

REFLECTIONS  FROM  SEA  ICE. 

During  the  past  seven  months,  Mr.  R.  C.  Byrd  has  measured  the 
complex  dielectric  constant  of  representative  samples  of  naturally- 
occurring  sea  ice,  as  a  function  of  temperature,  salinity,  and  orien¬ 
tation.  Details  are  reported  in  his  M. S.  thesis  of  December  1971, 
which  constitutes  the  appendix  of  this  report.  Preliminary  results 
were  presented  by  Mr.  Byrd  at  the  1971  IEEE  Conference  on  Engineering 
in  the  Ocean  Environment,  San  Diego,  September  21-24,  1971,  in  his 
paper  entitled,  "Interpretation  of  Radar  Returns  from  Sea  Ice." 

Sea  ice  may  be  considered  to  be  a  mixture  of  dielectrics,  if 
the  size  of  the  inhomogenieties  is  small  compared  to  a  wavelength. 

(Some  samples,  which  contained  entrapped  air  bubbles  of  dimensions 
comparable  to  a  quarter -wavelength,  demonstrated  erratic  behavior 
due  to  resonant-cell  scattering  effects.)  At  temperatures  wanner  than 
-21°C,  the  sea  ice  can  be  regarded  as  a  two-phase  mixture  of  relatively 
saline  brine  inclusions.  At  temperatures  colder  than  -21°C,  the 
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brine  inclusions  begin  to  freeze,  and  a  simplified  model  of  two 
solid  phases  (ice  crystals  with  frozen  brine  inclusions)  may  be 
used.  The  dielectric  loss  in  the  frequency  range  under  considera¬ 
tion,  26  -  40  GHz,  is  due  to  a  rotational  degree  of  freedom  of  the 
water  molecule  in  the  gaseous  or  liquid  state.  In  the  measurements 
reported  in  the  Appendix,  a  distinct  decrease  in  dielectric  loss 
was  observed  for  temperatures  colder  than  -21°C,  as  expected. 
Furthermore,  an  anistropy  was  observed  in  the  loss  tangent,  de¬ 
pending  upon  whether  the  electric  field  was  parallel  or  perpendi¬ 
cular  to  the  brine  inclusions.  This  anistropy  may  give  some  insight 
into  the  details  of  the  dielectric  loss  process.  However,  in  the 
case  of  a  radar  beam  striking  a  resolution  element  area  several  feet 
square,  a  random  combination  of  brine  inclusion  directions  at  the 
surface  is  likely,  so  that  the  polarization  anistropy  is  not  likely 
to  be  observed  in  the  radar  return  because  of  this  particular  effect. 
This  does  not  rule  out  the  possibility  of  an  anistropy  in  the  return 
due  to  other  causes;  future  field  measurements  will  determine  whether 
this  occurs. 

Some  insight  into  the  scattering  and  absorption  of  incident  radar 
waves  by  sea  ice  can  be  obtained  by  numerically  evaluating  the  atten¬ 
uation  distance  1/a  ,  the  distance  into  the  sea  ice  through  which  the 
wave  penetrates  before  it  decays  to  1/e  «  0.368  of  its  original  value. 
Von  Hippel1  has  given  the  expression  for  1/a  as: 

1  f  2 

a  "  2tt  •pTgflJTT  tan*  6  - 
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Assuming  an  average  frequency  of  34  GHz,  and  an  average  salinity 
of  4.48  °/oo,  the  dielectric  constant  ( z'/z is  3.08,  3.01,  and  3.30 
at  -7°C,  -25°C,  and  -32°C,  respectively,  while  the  loss  tangent  (tan  6) 
is  0.16,  0.028,  and  0.007,  respectively,  at  those  temperatures.  Num¬ 
erically  evaluating  Equation  1,  we  find  an  attenuation  distance  1/a 
of  1.03  cm.,  5.76  cm.,  and  22  cm.,  at  temperatures  of  -7°C,  -25°C, 
and  -32°C,  respectively.  A  common  winter  condition  is  a  temperature 
gradient  from  the  cold  ambient  air  temperature  at  the  surface  of  the 
sea  ice  to  the  -1.8°C  temperature  established  by  the  sea  at  the 
water-ice  boundary.  A  34  GHz  wavy  penetrating  the  ice,  launched  from 
an  airborne  SLAR  at  nearly  normal  Incidence,  would  only  penetrate  to 
a  depth  exceeding  a  meter  if  extensive  brine  drainage  had  occurred, 
and  the  air  temperature  was  very  cold.  Because  of  the  high  attenua¬ 
tion  of  sea  ice  at  temperatures  warmer  than  -10°C,  it  is  extremely 
unlikely  that  such  a  wave  could  penetrate  to  the  water-sea  ice  interface. 
An  abrupt  discontinuity  in  salinity  near  the  water-sea  ice  interface 
could  lead  to  a  reflected  wave,  however,  which  (if  detectable)  might 
give  a  measure  of  ice  thickness.  Further  analysis  and  field  measure¬ 
ments  are  required  to  resolve  this  question. 


for  the  electric  vector  is  (normal  incidence  case)1 
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where  a  wave  penetrates  from  a  material  of  intrinsic  impedance  (Z^) 
into  a  material  of  intrinsic  impedance  (Z^).  If  the  first  material 
is  free  space,  then  algebra  shows  that 


This  is,  in  general,  a  complex  number.  One  extreme  case  is  sea 

ice  that  is  several  years  old,  which  has  been  subjected  to  considerable 

brine  drainage  and  is  of  very  low  salinity.  Under  very  cold  winter 

conditions,  the  loss  tangent  of  the  ice  at  the  air -ice  Interface  may 

be  neglected.  Then  the  relative  dielectric  constant  is  near  3.1,  and 

|r  |  s  .275.  The  coefficient  of  energy  reflection  (R)  is 
£« 

R  -  |rE|2  (4) 

which  numerically  is  R  -  .0758,  for  the  particular  example  taken. 

This  represents  the  power  reflected  at  normal  incidence  from  very 

f 

old,  very  cold  ice.  A  pressure  ridge  in  winter  has  some  suall 
fractional  area  of  its  rough  surface  which  is  normal  to  the  incident 
wave,  and  the  backscatter  from  an  old  pressure  ridge  would  be  reduced 
by  this  effective  cross-section  factor.  This  model  of  backscattering 
is  drastically  oversimplified,  however;  when  one  encounters  higher 
temperatures,  the  loss  tangent  Increases,  increasing  R.  The  viewpoint 
of  "small  specular-reflecting  surface  area"  should  then  be  replaced 
by  a  more  complete  analysis  z* 3  which  shows  the  existence  of  a  diffuse 
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scattered  component  with,  in  general,  a  different  polarization  than 
the  incident  wave.  Field  measurements  will  be  taken  to  determine 
the  extent  of  this  effect,  and  correlation  of  that  data  with  Beckmann's 
theory  is  planned. 

It  appears  that  a  large  distinction  exists  between  an  open  or  a 
freshly  frozen  lead  and  older,  colder  flat  ice  surrounding  it,  when 
viewed  at  normal  incidence.  Future  analysis  will  quantify  this  dis¬ 
tinction. 

III.  FIELD  MEASUREMENT  PROGRAM 

As  indicated  in  the  preceding  section,  field  measurements  of 
backscatter  are  planned  for  sea  ice  of  varying  temperature,  salinity, 
and  topography.  Similar  measurements  are  planned  for  snow-covered  sea 
ice,  pure  ice,  and  soil,  so  that  a  comprehensive  picture  of  the  milli¬ 
meter-wave  radar  reflectivity  of  arctic  surfaces  may  be  known. 

During  the  past  seven  months,  several  representative  sites  for 
these  measurements  were  selected  and  photographed.  A  test  fixture 
was  fabricated  which  would  position  the  transmitting  and  receiving 
horn  antennas  exactly,  and  upon  which  they  could  be  moved  to  vary  the 
angle  of  incidence  of  the  wave. 

This  fixture  is  in  the  form  of  a  track  extending  over  a  90°  arc, 
upon  which  is  mounted  a  sliding  traveler  which  holds  the  antennas.  Both 
transmitting  and  receiving  horns  can  be  independently  rotated,  making 
possible  the  measurement  of  backscatter  for  horizontal,  vertical,  and 
crossed  polarization.  Considerable  ingenuity  was  embodied  in  the  fixture 
design  so  that  all  mechanical  positions  required  could  be  achieved 
easily  in  the  field,  with  positive  positioning.  No  mechanical 
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measurements  on  the  fixture  are  required  in  the  field,  reducing  sources 
of  error.  Tolerances  on  dimensions  and  angles  are  in  every  case  those 
obtainable  during  the  machining  process. 

A  simple  transmission  and  reflection  measurement  technique  is  used, 
in  which  an  r.f.  signal  is  amplitude-modulated  before  transmission,  and 
the  received  r.f.  energy  is  demodulated  and  detected  with  a  phase- 
locked  receiver.  A  received  signal  60  db  below  the  transmitted 
signal  could  easily  be  distinguished  using  this  approach,  which 
proved  adequate  to  detect  reflected  power  in  the  first  measurements 
on  snow-covered  soil. 

Complete  data  was  obtained  for  a  particular  site  near  the  end  of 
the  reporting  period.  This  site  had  approximately  28  inches  of  dry 
powder  snow  covering  the  soil,  with  an  ambient  temperature  in  the  range 
of  -25°C  to  -45°C.  Analysis  of  the  data  is  in  progress,  and  the  data 
will  be  presented  and  discussed  fully  in  the  next  technical  report. 

During  the  next  quarter,  data  will  be  obtained  on  sea  ice  at 
Barrow,  and  on  snow-covered  ice  and  terrain  at  other  locations  near 
Barrow  and  Fairbanks. 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  dielectric  properties  of  sea  ice  have  been  measured  in  the 
frequency  range  26-40GHz.  A  comprehensive  report  on  this  subject  is 
included  as  the  Appendix.  The  fixtures  for  field  measurements  have 
been  designed,  fabricated,  and  tested.  The  backscatter  from  snow- 
covered  soil  at  one  site  has  been  measured  as  a  function  of  incidence 
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angle  and  polarization. 

During  the  next  reporting  period,  additional  field  measurements 
of  backscatter  will  be  taken  for  sea  ice,  fresh-water  ice  and  land, 
both  snow-covered  and  exposed.  Laboratory  measurements  of  the 
dielectric  properties  of  snow  and  pure  ice  will  commence.  A  contin¬ 
uation  of  this  research  is  recommended  to  measure  backscatter  during  the 
summer  and  early  winter,  and  to  analyze  field  measurements.  This 
should  lead  to  an  understanding  of  the  relative  contributions  of  the 
surface  roughness  and  of  the  surface  dielectric  properties. 
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ABSTRACT 


The  dielectric  properties  of  sea  ice  are  measured  for  fre¬ 
quencies  from  26  GHz  to  40  GHz.  The  effect  of  temperature  and 
salinity  on  the  complex  permittivity  are  determined.  The  orienta¬ 
tion  of  the  brine  cells  within  the  sea  ice  is  found  to  have  an 
effect  on  the  loss  tangent  of  the  ice.  An  average  value  for  the 
real  part  of  the  complex  permittivity  is  derived  along  with  an 
empirical  expression  for  determining  loss  tangent. 


v 


I.  INTRODUCTION 


With  the  Increase  In  interest  in  recent  years  in  recovering 
natural  resources  from  the  Arctic  region,  new  and  better  methods 
of  transport  through  and  over  the  Arctic  Ocean  ice  cover  are  being 
sought.  To  facilitate  routing  of  icebreakers  and  surface-effect 
vehicles,  better  ice  navigation  systems  and  methods  of  predicting 
ice  surface  topography  are  needed.  One  area  of  major  emphasis  is 
in  the  development  of  radar  systems  which  produce  a  graphic  repre¬ 
sentation  of  the  ice  surface,  both  for  large  area  prediction,  as 
in  the  case  of  airborne  side-looking  radars,  and  for  navigational 
radar  carried  by  the  vessels  themselves.  Because  of  the  nature 
of  the  surface  to  be  represented  and  the  information  that  is 
being  sought  from  the  radar  presentations,  high  degrees  of  resolu¬ 
tion  are  desired.  This  will  necessitate  use  of  systems  at  fre¬ 
quencies  in  excess  of  20  GHz.  One  problem  which  has  not  been 
satisfactorily  solved,  involving  the  design  of  sea  ice  radar  systems 
is  development  of  a  method  of  interpreting  the  return  signals 
received  by  such  systems.  A  radar  Impulse  returning  from  a  sea 
surface  having  partial  ice  cover  may  represent  reflections  from 
fresh  water,  fresh  water  ice,  sea  water,  and  sea  water  ice  in  a 
variety  of  combinations — all  having  somewhat  different  reflectivity 
characteristics.  The  return  slgual  is  the  result  of  a  complex 
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interaction  of  the  electromagnetic  field  of  the  impulse  with  the  sea 
ice  profile.  The  dielectric  constant  and  loss  tangent  of  the  sea 
cover  components  partly  determine  the  reflection  coefficient  at 
each  boundary.  The  differences  in  reflective  characteristics  be¬ 
tween  different  sea  ice  profiles  will  be  better  understood  when 
more  information  is  available  about  the  dielectric  properties  of 
sea  ice.  This  study  seeks  to  determine  these  properties  for  natural 
sea  ice  in  the  frequt  _y  from  26.5  GHz  to  40.0  GHz,  and  their  changes 
with  salinity  and  temperature. 
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II.  GENERAL  PROPERTIES  OF  SEA  ICE 


In  this  section  the  physical  properties  of  sea  ice  which 
effect  its  electromagnetic  characteristics  are  discussed,  along 
with  a  brief  description  of  the  relationship  between  reflectivity 
and  the  dielectric  properties  of  the  reflecting  material. 

1.  Physical  Properties  of  Sea  Ice 

Sea  ice  differs  from  fresh  water  ice  in  that  it  has  liquid 
brine  inclusions  contained  in  its  ice  matrix.  The  rate  at  which 
the  freezing  process  takes  place  determines  to  a  great  extent  the 
amount  of  salts  trapped  in  the  ice,  and  therefore,  the  percent  of 
brine  volume  to  be  found.1  Rapid  freezing  produces  more  included 
brine.  As  the  temperature  continues  to  decrease,  the  relative 
volume  of  the  included  brine  decreases  as  more  water  freezes  out, 
with  an  increase  in  salt  concentration  in  that  brine  which  remains. 
At  about  -22°C  solid  salts,  predominantly  NaCl,  begin  to  precipitate 
out  of  the  brine.  2  An  example  of  the  formation  pattern  of  sea  ice 
is  shown  in  Figure  1,  illustrating  the  approximate  arrangement  of 
the  brine  inclusions  in  the  ice  matrix.  The  model  was  used  by 
Assur  and  Weeks3  to  simulate  the  mechanical  properties  of  sea  ice 
and  is  greatly  oversimplified  with  regard  to  the  orderly  arrangement 
of  the  brine  inclusions  and  their  shape.  The  inclusions  have  been 
observed  to  sometimes  be  connected  in  a  direction  perpendicular 
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to  their  vertical  axis  and  their  cross-sectional  shape  is  generally 
not  as  well  defined  as  that  shown. 

Once  the  ice  has  reached  its  minimum  temperature,  the  process 
has  by  no  means  stabilized.  The  liquid  brine  inclusions  tend  to 
gravity  drain  to  a  lower  point  and  eventually  out  of  the  ice.  If 
the  ice  cycles  through  warming  periods,  the  amount  of  trapped  brine 
decreases  more  rapidly  due  to  increased  migration. **  The  maximum 
salinity  found  in  young  sea  ice  is  generally  around  20°/oo  5  (parts 
per  thousand)  and  it  may  range  to  nearly  zero  in  old  ice.  The 
average  salinity  in  the  sea  ice  which  comprises  he  majority  of 
ice  surface  disruptions  is  around  5°/oo.6 

In  predicting  the  electromagnetic  properties  of  sea  ice  at 
microwave  frequencies,  the  properties  of  its  two  major  constituents, 
ice  and  liquid  brine,  must  be  considered.  Pure  ice  has  been  found 
to  be  a  low- loss  dielectric  at  frequencies  above  107  Hz,  while  water, 
which  the  brine  is  expected  to  resemble,  has  been  predicted  to  have 
high  loss  in  the  frequency  range  between  10* 0  and  1011  Hz.7  The 
percent  of  brine  volume  in  sea  ice  varies  greatly  with  salinity  and 
temperature;  Figure  2  shows  this  relationship.  2  The  dielectric 
properties  of  a  particular  sample  of  sea  ice  will  depend  on  its 
past  history  and,  more  important,  on  its  current  temperature. 

2.  Electromagnetic  Properties  of  Dielectric  Materials. 

Before  proceeding  with  the  discussion  of  the  measurement  of 
the  electromagnetic  properties,  we  shall  examine  briefly  the  general 
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V 

properties  of  dielectric  materials  and  how  they  relate  to  the  re¬ 
flectivity  of  electromagnetic  waves. 

A  simple  electromagnetic  wave  traveling  in  free  space  can  be 
visualized  as  consisting  of  two  components,  as  shown  in  Figure  3; 
an  electric  field  vector  (E)  and  a  magnetic  field  vector  (H).  In 
the  undisturbed  case  these  vectors  are  at  right  angles  to  each 
other  and  each  is  at  right  angles  to  the  direction  of  propagation. 

Our  problem  is  to  define  the  information  necessary  to  determine  how 
the  wave  reacts  when  confronted  with  substances  other  than  free  space. 
When  this  occurs,  the  incident  wave  is  broken  into  two  additional 
waves;  the  reflected  wave — that  portion  which  is  returned  from  the 
boundary  of  the  dielectric,  and  the  transmitted  wave — that  portion 
which  proceeds  through  the  dielectric.  The  three  waves  are  inter¬ 
connected  by  the  boundary  condition  that  the  tangential  component 
of  E  and  H  must  be  continuous  in  traversing  the  Interface  between 
two  media.  This,  in  fact,  is  two  conditions,  since  the  amplitudec 
as  well  as  the  phases  must  be  continuous.  Figure  4  shows  the  simple 
case  of  a  plane  wave  at  normal  incidence  from  one  dielectric  media 
to  another.  We  are  generally  not  interested  in  the  absolute  value 
of  the  reflected  or  transmitted  wave  amplitudes,  but  are  interested 
in  the  ratios  of  the  different  amplitude  components,  the  reflection 
coefficient  (r)  and  the  transmission  coefficient  (t).  These  may  be 
defined  in  terms  of  the  intrinsic  impedances  of  the  media.6 


The  intrinsic  impedances  may  in  turn  be  defined  in  terms  of  the 
complex  permittivity  (e*)  and  the  complex  permeability  (y*),8 


where  the  complex  propagation  constant  (y)  is  defined 

y  =  jai(e*y*)J4  -  a  +  jB  (4) 

where:  to  =  radian  frequency  of  the  wave. 

Figure  5  illustrates  a  wave  at  other  than  normal  incidence. 

The  relationship  of  the  different  wave  components  to  intrinsic 
impedances  and  angles  of  incidence  and  refractions  are  given  by 
von  Hippel.8  (see  the  equations  contained  in  drawing  #5  as 
illustration).  These  figures  and  expressions  are  somewhat  simplified, 
but  serve  to  illustrate  the  point  that  the  characteristics  of  electro¬ 
magnetic  waves  reflected  and  refracted  by  dielectrics  are  dependent 
on  the  complex  permittivity  and  permeability  of  the  substances  on 
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which  the  waves  are  incident.  Most  dielectrics  o£  interest  are 

for  all  practical  purposes  non-magnetic,  therefore,  y*  *  p  ,  the 

o 

permeability  of  free  space.  Therefore,  the  first  step  in  under¬ 
standing  the  reflective  characteristics  of  a  particular  dielectric 
material  comes  with  knowing  itd  complex  permittivity  (e*)  for  the 
frequencies  of  interest. 


-12- 


III.  COMPLEX  PERMITTIVITY  DETERMINATION 


Methods  for  measuring  the  dielectric  properties  of  materials 
at  high  frequencies  have  been  described  by  a  number  of  authors  and 
vary  greatly  in  the  specific  details,  depending  on  the  range  of 
accuracy  for  which  they  may  be  used  and  the  physical  properties 
of  the  materials  to  be  measured.8'  9#  10/  11  Most  of  the  methods 
described  are  intended  for  use  in  homogeneous  materials  where  a 
high  degree  of  accuracy  is  possible;  they  require  a  high  degree  of 
control  of  sample  dimensions,  measurement  of  absolute  distances,  and 
generally  require  complicated  calculations  for  each  measurement 
taken.  A  method  which  has  been  found  to  be  particularly  well  suited 
for  measurements  in  sea  ice  at  the  frequencies  to  be  considered  here 
was  described  by  Tinga  and  Edwards;  a  swept  frequency  technique, 
it  offers  several  advantages  over  other  methods  available.11  This 
technique  requires  only  one  length  measurement,  that  of  the  sample. 

It  provides  for  rapid  acquisition  of  data  and  can  be  used  in  samples 
having  some  inhomogeneities.  It  offers  the  opportunity  for  a  trade 
of  some  accuracy  for  a  great  deal  of  simplicity. 

Most  conventional  methods  of  dielectric  measurements  are  sim¬ 
plified  considerably  if  the  sample  length  has  a  particular  value 
such  as  an  odd  or  even  number  of  quarter  wave-lengths  in  the  material. 
However,  since  wavelength  in  the  material  is  a  function  of  the 
dielectric  constant,  this  condition  is  usually  difficult  to  achieve 
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in  practice,  since  this  is  the  characteristic  to  be  measured.  The 
electrical  length  of  the  sample,  however,  can  be  varied  easily  by 
changing  the  frequency.  If  the  assumption  is  made  that  the  dielec¬ 
tric  properties  of  the  material  under  test  do  not  change  greatly 
with  frequency  over  a  relatively  short  frequency  range,  then 
changing  the  electrical  length,  by  sweeping  the  frequency  over 
some  band,  has  the  effect  of  changing  the  physical  length  of  the 
sample.  This  allows  a  situation  to  be  created  which  makes  it 
relatively  easy  to  calculate  the  complex  permittivity  of  the 
material.  This  approach  forms  the  basis  for  the  swept  frequency 
technique.  The  experimental  measurement  is  basically  a  reflection 
measurement,  which  relies  on  the  interface  reflections  and  the 
attenuation  of  the  sample  to  give  information  about  the  complex 
permittivity. 1 1 

A  block  diagram  of  the  experimental  apparatus  is  shown  in 
Figure  6.  A  leveled  microwave  sweep  generator  sends  an  incident 
wave  to  the  sample,  which  will  reflect  and  transmit  some  of  the 
energy.  If  the  total  reflected  signal  is  detected  and  displayed 
versus  frequency,  a  characteristic  pattern  such  as  is  given  in 
Figure  7  results.  This  pattern  contains  a  large  amount  of  infor¬ 
mation  about  the  dielectric  behavior  of  the  material  under  exam¬ 
ination. 

The  following  is  an  outline  of  the  derivation  of  the  necessary 
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Reflection  Ratio 


x-y 


Plotter 


Fig.  6  Equipment  Setup 


Fig.  7  Representative  Data  Trace 


equations,  taken  from  the  detailed  work  of  Tinga  and  Edwards,  for 
the  particular  case  of  a  short  circuit  sample  holder  termination, 
which  was  used  in  these  measurements  and  is  shown  in  Figure  8. 

The  propagation  constant  for  a  medium  of  arbitrary  dielectric 
constant,  e*,  filling  a  rectangular  wave-guide,  may  be  expressed 


y  =  a  +  jg  =  J(2tt  f/c)  ^  e*-  P 

where 


(7) 


e*  =  t'  +  ]e"  (complex  permittivity)  (8) 

P  *  (f  /f)2  (9) 

c 

with  f  being  the  cutoff  frequency  of  the  wave-guide  used  and  f  as 
c 

the  frequency  being  examined;  c  represents  the  speed  of  light  in 
free  space. 

Squaring  (7)  and  solving  for  t'  and  z"  gives 

e'  =  P  +  (c/2irf ) 2  (  a2-  g2)  (10) 

e"  -  2(c/2TTf)2  ag  (11) 


To  find  values  for  a  and  g,  referring  to  Figure  8, 


R  =  -  = 
A 


rl  +  r^e 


— 2y 


1  +  rir2e 


rz 


(12) 


where 

A  =  amplitude  of  the  incident  wave  electric  field 

B  =  amplitude  of  the  total  reflected  wave  electric  field 

rx  =  reflectivity  coefficient  at  face  of  the  sample 

r2  =  reflectivity  coefficient  at  the  short  circuit 

l  =  length  of  the  sample 

Yj  =  propagationconstant  of  the  sample 
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For  the  case  of  a  short  circuit  termination,  r  =  -1,  and 

2 

r  =  /  1  -  P  -  ^e*-  P 

1  - - - 

-  P  +  ~  P  (13) 

To  find  the  points  of  extrema  in  equation  (12),  the  condition 
that  the  first  derivative  with  respect  to  frequency  be  equal  to  zero 
is  applied.  To  find  a  first  approximation  to  the  solution,  the  assump¬ 
tions  of  plane  wave  propagation  and  a  lossless  dielectric  (  a=  0)  are 
made,  resulting  in 


B  -  Yi  n  =  integer,  for  extrema  (14) 

A3  =  m  =  integer  (15) 

Assuming  that  a2<<32  one  may  solve  equation  (10)  for  3  • 

3  =  (2tt  f/c)Ve'  -  P  (16) 

* 

Assuming  that  e  remains  constant  between  two  adjacent  maxima  or 
minima  and  taking  m  =  2,  since  the  electrical  length  between  two  ad¬ 
jacent  maxima  or  minima  equals  one-half  the  wave  length  within  the 
sample,  substitituting  equation  (16)  into  equation  (15)  yields 

(c/2«,Af)  =  z'l/P~  (17) 

1  1 
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where  Af  =  f  -  f  (See  Figure  7) 
2  1 

Squaring  and  solving  for  , 


£■'  =  1/2  E  [1  +  /( 1  -  4P)/E  ]  (18) 

where 

P  =  (f2p2  +  fiPi) 

fi  +  f2 

E  =  (c/2£Af )  2 


It  can  be  seen  in  equation  (18)  that  all  variable  terms  are  available 
from  the  experimental  trace  and  that  c'  is  a  function  of  Af  only. 

Now,  taking  equation  (12)  and  assuming  that  r  t  has  only  a  real 


component;  i.e.,  the  sample  is  lossless, 

a  .  b-l- 


1  -  r,  e 


-20,1 


(19) 


Solving  equatic..  (19)  for  a  gives 

a  =  1/2S,  (InA)  (20) 

where 

A  =  (1  -  R|rJ)(R  -  |rj)  (21) 

The  bars  over  the  ratios  indicate  that  they  are  taken  for  maximum  values. 

Substituting  equations  (20)  and  (16)  into  equation  (11) 

e"  =  (c/2TTf  i)  (  J  e'  -  P  )  In  A  (22) 

1  m  '  1 

or  _ 

tan  6  =  (c/ 2tt f  e'  )  (Jz'  -  P)  In  A  (23) 

m  1  1 

Equations  (18)  and  (23)  represent  the  information  necessary  to  deter¬ 
mine  the  complex  permittivity  for  materials  which  approximately  obey  the 
assumptions  made  in  arriving  at  these  solutions. 

In  the  formation  of  the  expression  for  the  dielectric  constant  (c') , 
equation  (18),  plane  wave  propagation  along  with  a  lossless  dielectric 
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was  assumed.  These  conditions  are  not  generally  the  case  in  practical 
laboratory  experiments.  However,  Tinga  and  Edwards  have  shown  by  use 
of  computer  solutions  for  wave-guide  propagation  in  lossy  dielectrics 
that  these  result  only  in  a  phase  shift  in  the  value  of  $  expressed  in 
equation  (14),  but  do  not  affect  the  difference  expression  as  stated  in 
equation  (15),  for  a  given  sample  length.  The  major  assumption  left 
to  be  fulfilled  in  satisfying  equation  (18)  is  that  the  dielectric 
constant,  e' ,  not  vary  greatly  between  any  two  minima,  since  this 
expression  gives  the  average  value  of  this  property  between  the  frequen¬ 
cies  where  Af  is  measured.  The  results  to  be  presented  later  will  show 
that  this  is  a  valid  assumption  for  sea  ice  in  the  frequency  range  con¬ 
sidered  and  that,  in  fact,  the  dielectric  constant  generally  remains 
within  -  10  percent  of  the  same  value  over  the  entire  band  measured, 
for  most  conditions. 

In  arriving  at  the  expression  for  tan  6,  equation  (23),  it  was 
assumed  that  a2<<02,  which  is  a  good  assumption  for  dielectric  mater¬ 
ials  to  which  this  method  can  be  applied,  since  a  trace  as  shown  in 
Figure  7  cannot  be  achieved  for  real  sample  lengths  unless  this  con¬ 
dition  is  met.  It  was  also  assumed  that  r  could  be  found  using  only  the 
real  part  of  the  complex  permittivity  in  arriving  at  a  value  for  A  in 
equation  (21).  This  will  result  in  some  error  if  the  value  of  tan  6 
is  large,  requiring  reinteration  using  the  full  value  of  permittivity 
for  greater  accuracy.  However,  for  values  of  tan  6<0.1,  the  percentage 
error  of  the  calculations  using  only  the  real  part  of  the  complex 
permittivity  in  calculating  r  has  been  found  to  be  less  than  5%. 11 
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This  is  generally  the  case  for  sea  ice  where  a  usable  sample  trace  is 
achieved.  It  was  found  that  when  the  value  of  loss  tangent  was  much 
above  0.1,  the  minima  could  not  be  defined  on  the  sample  trace.  This 
occurs  only  at  high  salinities  and/or  high  temperatures,  as  will  be 
shown,  and  an  alternate  method  must  be  used  for  future  laboratory 
measurements  of  very  lossy  sea  ice. 
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IV.  EXPERIMENTAL  PROCEDURE 


Ice  samples  for  this  experiment  were  taken  from  the  Arctic  Ocean 
in  January  and  March  of  1971  from  first-year,  shore-fast  ice.  Care 
was  taken  to  insure  that  the  ice  was  kept  at  temperatures  below  -40°C 
until  the  time  of  use.  The  ice  was  selected  at  random  from  undisturbed 
pieces  of  pack  ice  at  locations  ranging  from  on  the  surface  (under  snow 
cover)  to  eighteen  inches  beneath  the  surface.  All  samples  fell  into 
three  groups  by  salinity:  2.85  °/oo,  3.4  °/oo,  and  7.2  °/oo. 

A  3db  directional  coupler  was  selected  as  the  ice  sample  holder 
illustrated  in  Figure  9.  As  indicated,  the  straight,  uncoupled  end 
section  of  the  coupler  was  used  as  the  holder.  This  provided  a  device 
in  which  the  wave  propagation  was  understood  and  allowed  the  number 
of  components  used  in  the  experiment  to  be  kept  at  a  minimum.  A 
copper  shorting  plate  was  used  as  a  termination. 

A  number  of  temperature  control  mechanisms  were  considered  and 
the  following  arrangement  was  found  to  be  most  convenient: 

1)  A  temperature-sensing  thermistor  was  mounted  at 
the  approximate  mid-point  of  the  sample  on  the 
sample  holder. 

2)  The  thermal  mass  of  the  sample  holder  was  increased 

by  addition  of  a  tightly-wound  casing  of  fine  copper  wire. 

3)  The  entire  apparatus  was  covered  with  polyurethane  insul¬ 
ation  approximately  6  inches  thick,  so  that  the  temperature 
decay  from  -60°C  to  0°C  when  exposed  to  ambient  room  tempera¬ 
ture  averaged  less  than  1.5°C  per  minute,  and  was  considerably 
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3)  continued... 

slower  at  the  warmer  temperatures  where  the  dielectric 
properties  were  expected  to  change  more  rapidly. 

Before  each  series  of  measurements,  the  sample  and  sample  holder 
were  allowed  to  stabilize  at  -60°C.  The  sample  holder  was  then  re¬ 
moved  from  cold  storage  and  allowed  to  slowly  warm,  with  readings 
being  taken  at  intervals  between  -45°C  and  -7°C.  The  sample  holder 
was  returned  to  cold  storage  without  allowing  the  sample  to  melt. 

A  swept  backward  wave  oscillator,  Hewlett-Packard  Model  8690B, 
was  used  as  the  signal  generator.  A  power  level  of  approximately 
5  milliwatts  incident  upon  the  sample  was  used  for  all  cases.  An 
automatic  leveling  control  loop  was  used  in  the  sweep  oscillator 
to  insure  that  the  magnitude  of  the  output  power  was  the  same  at 
all  frequencies.  The  band  from  26.5  GHz  to  40.0  GHz  was  swept  in 
approximately  five  seconds  on  all  samples.  This  sweep  time  was 
found  to  be  slow  enough  for  accurate  following  by  the  X-Y  plotter, 
Hewlett-Packard  Model  7004B,  but  fast  enough  so  that  less  than  0.2°C 
change  in  temperature  occurred  during  the  measurement. 

The  reference  for  determining  the  value  of  the  incident  wave 
electric  field  was  taken  from  the  reflected  wave  of  the  shorted 
sample  holder  without  a  sample  in  place.  This  allowed  for  cancella¬ 
tion  of  effects  due  to  the  directional  coupler.  Reference  traces 
were  run  immediately  before  the  sample  holder  was  chilled  and  again 
after  the  sample  was  removed,  to  check  drift  of  the  oscillator. 
Reference  traces  were  also  run  on  the  empty  sample  holder  while  warming 
from  -60°C  to  determine  the  effect  of  temperature  change;  none  was 
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noted.  However,  it  was  determined  that  the  crystal  detectors  used 
to  detect  the  magnitude  of  the  reflected  power  were  sensitive  to 
temperature  and  had  to  be  protucted  from  temperature  change  in  order 
to  insure  consistency.  This  was  done  by  allowing  the  detectors  to 
come  in  contact  with  the  cold  apparatus  only  during  the  short  inter¬ 
vals  necessary  to  take  readings.  Removing  and  attaching  the  detectors 
had  no  noticeable  effect  on  the  measurements. 

The  individual  sample  data  was  recorded  on  an  X-Y  plotcer,  as  shown 
in  Figure  6.  Figure  10  shows  a  data  trace  for  a  typical  sample. 

Usable  data  traces  were  received  from  most  samples  at  temperatures 
of  -32°C,  -25°C,  -21.5°C,  -16.5°C,  and  -7°C.  Below  -32°C  and  above 
-7°C  the  loss  factor  of  the  sea  ice  examined  was  not  within  the  range 
of  this  detection  technique. 

The  most  difficult  and  time-consuming  portion  of  this  experiment 
was  in  preparing  the  sea  ice  samples  for  the  sample  holder.  As  stated 
earlier,  care  was  taken  to  prevent  the  samples  from  being  exposed  to 
temperatures  above  -40°C  for  any  appreciable  length  of  time  in  order  to 
prevent  brine  drainage.  The  raw  sea  ice  and  prepared  samples  were 
stored  at  -60°C  in  the  laboratory  while  awaiting  use.  Prepared 
samples  were  sealed  in  airtight  containers  to  prevent  evaporation 
in  storage.  Although  a  number  of  methods  of  actually  preparing  samples 
were  tried,  no  efficient  method  was  found,  and  in  the  end,  samples 
were  prepared  by  manually  carving  with  a  razor  blade.  A  simple  jig  was 
used  to  insure  square  corners,  and  a  proper  fit  in  the  wave-guide.  The 
samples  were  rectangular  in  cross-section,  to  match  the  wave-guide,  and 


* 
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ranged  in  length  from  1.2  to  3.5  cm.  Care  was  taken  in  selecting 
the  samples  so  that  they  contained  as  few  cracks,  bubbles  and  similar 
inhomogeneities  as  possible.  The  samples  were  examined  under  a  micro¬ 
scope  during  preparation  so  that  they  could  be  oriented  with  the 
brine  inclusions  either  perpendicular  or  parallel  to  the  direction  of 
the  electric  field  in  the  wave-guide.  Only  samples  with  well-defined 
and  reasonably  consistent  brine  inclusions  were  used.  Tb  samples 
were  prepared  in  a  cold  room  at  -22°C.  Care  was  taken  to  minimize 
the  time  that  the  samples  were  subjected  to  this  temperature  to 
prevent  brine  drainage.  Sample  preparation  time  averaged  about 
1-1/2  hours  for  each  sample.  Thirty  samples  were  prepared,  twenty 
were  successfully  examined. 


V.  RESULTS 

Tables  AI,  All,  and  AIII  of  the  Appendix  show  the  values  of  indi¬ 
vidual  data  points  obtained  in  this  experiment  for  the  three  sea  ice 
salinity  groupings  respectively.  Figures  Al  through  A30  present  a 
least  squares,  best  straight  line  fit  of  the  tabulated  data  points  for 
the  different  salinities,  temperatures  and  sample  orientations  examined 
and  an  average  value  for  the  individual  conditions  without  regard  to 
sample  orientation.  It  will  be  noted  that  in  some  cases  data  is  not 
available  for  parallel  sample  orientation  at  -7°C  and  perpendicular 
sample  orientation  at  -32°C.  Most  of  the  samples  measured  were 
too  lossy  to  be  measured  in  the  former  condition  and  not  lossy  enough 


for  the  latter  condition.  The  range  over  which  samples  could  be 
evaluated  extended  from  a  loss  tangent  of  approximately  0.120  down 
to  0.001. 

In  evaluation  of  the  tabulated  data  points,  it  should  be  kept  in  mind 
that  each  point  represents  an  average  value  of  the  complex  permittivity 
over  a  small  band  of  frequencies,  generally  about  3  x  10  9  Hz,  and  that 
the  frequency  value  listed  represents  the  mid-point  of  that  band;  for 
example,  a  point  listed  f^.  28  x  109  represents  the  average  value 
from  26.5  to  29.5  x  109  Hz. 

Tables  I  through  V  show  the  values  of  the  real  part  of  the  dielec¬ 
tric  constant  (£.')  taken  from  the  average  value  plots  in  the  Appendix 
with  respect  to  temperature  and  salinity  for  even  frequencies  across 
the  band  measured.  In  Table  I  the  values  listed  for  the  7.20  °/oo 
salinity  samples  are  those  for  the  perpendicular  orientation  of  the 
brine  inclusions  only.  The  parallel  oriented  samples  were  too  lossy 
to  be  measured  fs  noted  above.  The  assumption  is  made  that  the  effect 
of  orientation  on  the  real  part  of  the  complex  permittivity  is  small. 

In  Table  V,  the  values  listed  for  the  2.85  °/oo  salinity  samples  are 
for  parallel  orientation  only.  Again,  the  assumption  of  a  small  change 
iue  to  orientation  is  made.  All  other  values  listed  are  for  the 
average  of  the  two  orientations  measured  at  a  particular  salinity, 
temperature,  and  frequency,  taken  from  the  data  plots  in  the  Appendix. 
Tables  VI  shows  the  average  values  of  the  real  part  of  the  dielectric 
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constant  with  respect  to  frequency  and  salinity  alone.  The  results 
of  these  tables  are  presented  graphically  in  Figures  11,  12,  and  13. 

Tables  VII  through  XI  show  the  values  of  loss  tangent  (tan  6) 
taken  from  the  average  value  plots  in  the  Appendix  with  respect  to 
temperature,  salinity  and  orientation  for  the  even  frequencies.  In 
table  VII  the  values  listed  for  the  average  value  of  loss  tangent  at 
-7°C  where  estimated  from  the  observation  that  at  the  warmer  tempera¬ 
tures,  the  average  loss  approaches  twice  the  value  of  that  of  the 
perpendicular  orientation.  The  values  listed  for  the  parall^ 
orientation  were  estimated  from  a  four  to  one  ratio  between  parallel 
and  perpendicular  orientation.  Examination  of  the  plotted  data  in  the 
Appendix  will  show  that  this  appears  to  be  a  reasonable  assumption. 
Table  XII  shows  the  average  loss  tangent  with  respect  to  frequency  for 
salinity  and  orientation.  These  tables  are  presented  graphically  in 
Figures  14,  15,  16,  and  17. 

The  data  which  was  used  to  produce  these  graphs  was  subject  to 
error  which  should  be  considered  in  their  interpretation.  Error  due 
to  the  assumption  of  constant  real  part  of  the  complex  permittivity 
with  frequency  made  in  arriving  at  equation  (18)  is  calculated  to  be 
less  than  1%  for  a  given  data  point.  The  assumption  of  low  loss  made 
in  arriving  at  equation  (23)  causes  a  5 %  error  in  loss  tangent  at 
-7°C  and  less  than  1%  at  -32°C.  This  error  will  cause  the  data  to 
be  lower  than  the  actual  value.  The  greatest  source  of  error  comes 
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from  interpreting  the  data  traces.  The  accuracy  varies  with  temperature 
from  -  5%  for  -7°C  through  -25°C  and  -  40%  for  values  at  -32°C.  These 
errors  are  random  in  nature  and  tend  to  be  cancelled  by  the  averaging 
process  used  to  arrive  at  the  straight  line  plots  shown  in  the  Appen¬ 
dix,  however,  they  are  apparent  in  examination  of  the  tabulated  data 
points  shown  in  Tables  AI  through  AIII  of  the  Appendix.  Error  due  to 

the  accuracy  of  the  sample  length  measurement  is  -  1%  for  the  real  part 

+ 

and  -  2%  for  the  loss  tangent. 

Overall  accuracy  is  expected  to  be  approximately  -  2%  for  the  real 
part  of  the  complex  permittivity  and  -  5%  for  the  loss  tangent. 
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TABLE  I  Comparison  of  the  Average  Values 
of  the  Real  Part  of  the  Complex  Permittivity 


Temperature:  -7°C 

9 

Frequency  (xlO  hz) 

;r  ro  32  34  36  38  40  Average 


3.00  3.07  3.13  3.19  3.26  3.32  3.39  3.19 


3.00  3.03  3.07  3.10  3.13  3.17  3.20  3.10 


3.23*  3.14*  3.05*  2.95*  2.85*  2.75*  2.65*  2.95* 


4J 

5  3.40 


nj 

w  7.20 


Average  3.08  3.08  3.08  3.08  3.08  3.08  3.08  3.08 


*values  for  perpendicular  orientation  only 


TABLE  II 


Temperature:  -16.5°C 


Frequency  (xlO^Hz) 

28  30  32  34  36  38  40  Average 


•5  3.40 

r-i  —  » ■ 

CO 

W  7.20 


2.99 

3.00 

3.01 

3.02 

3.03 

2.88 

2.93 

2.98 

3.03 

3.08 

3.04 

3.02 

2.99 

2.97 

2.95 

Average  2.97  2.98  2.99  3.01  3.02  3.03  3.04  3.01 
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Comparison  of  the  Average  Values  of  the  Real 
Part  of  the  Complex  Permittivity  (Continued) 
TABLE  III 


Temperature:  -21.5°C 

Frequency  (xlO^Hz) 

28  30  32  34  36  38  40  Average 


Salinity 

2.85 

2.86 

2.92 

2.98 

3.04 

3.10 

3.15 

3.20 

3.03 

3.40 

2.91 

2.93 

2.96 

2.98 

3.00 

3.03 

3.05 

2.98 

7.20 

2.86 

2.88 

2.90 

2.93 

2.96 

2.98 

3.00 

2.93 

Average 

2.88 

2.91 

2.95 

2.98 

3.02 

3.05 

3.08 

2.98 

TABLE  IV 


Temperature:  -25°C 

Q 

Frequency  (xlO  Hz) 


28 

30 

32 

34 

36 

38 

40 

Average 

>* 

2.85 

2.88 

2.93 

2.98 

3.03 

3.08 

3.13 

3.18 

3.03 

4J 

•H 

C 

•H 

3.40 

2.98 

2.98 

2.97 

2.97 

2.97 

2.96 

2.96 

2.97 

•H 

to 

7.20 

2.83 

2.90 

2.96 

3.02 

3.08 

3.14 

3.20 

3.02 

Average 

2.90 

2.94 

tm 

3.01 

3.04 

3.08 

3.11 

3.01 
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Comparison  of  the  Average  Values  of  the  Real 
Part  of  the  Complex  Permittivity  (Continued) 


TABLE  V 

Temperature:  -32°C 

Frequency  (xlO^Hz) 


28 

30 

32 

34 

36 

38 

40 

Average 

2.85 

3.25*  3.21* 

3.16* 

3.11* 

3.06* 

3.01* 

2.96* 

3.11* 

3.40 

2.73 

2.82 

2.91 

3.00 

3.08 

3.17 

3.25 

2.99 

7.20 

2.58 

2.99 

3.40 

3.80 

4.20 

4.59 

4.98 

3.79 

Average  j|  2.85  3.01  3.16  |  3.30  |  3.45  3.59  | 3.73  |  3.30 

*values  for  parallel  orientation  only 

TABLE  VI 

c'  versus  Frequency  (x109Hz) 


Average 

2. 85%o 

3.40%o 

7 . 20%o 

f 

e' 

> 

e 

£■' 

f' 

28 

2.94 

3.00 

2.90 

2.91 

30 

2.98 

3.03 

2.94 

2.99 

32 

3.03 

3.05 

2.98 

3.06 

34 

3.08 

3.08 

3.02 

3.13 

36 

3.12 

3.11 

3.06 

3.21 

38 

3.17 

3.13 

3.09 

3.28 

40 


3.21 


3.16 


3.12 


3.34 


Comparison  of  the  Average  Values 
of  Loss  Tangent 


TABLE  VII 


Temperature:  -7°C 


Frequency  (xlO^Hz) 


28 

30 

32 

34 

36 

38 

40 

Average 

m 

.268* 

.280* 

.292* 

.300* 

.308* 

.320* 

.328* 

.300* 

2.85 

AV. 

HH 

mm 

.150* 

0 

L 

■EM 

TflrTil 

WiiYkkk 

mi 

KEH 

mil 

IHQKEI 

wmest a 

.308* 

.320* 

.332* 

.348* 

wtia»Ei 

u 

3.40 

AV. 

■m 

— 

HI 

mm 

.  wm 

■B 

•H 

G 

L 

.070 

mm 

B-*  ■  i 

mi 

.080 

•H 

rH 

umn 

IfcIHiLl 

■Bra 

PIcfcMJ 

wmwK 

.360* 

.372* 

BOB 

in 

7.20 

AV. 

1 

■■ 

.172* 

L 

\mtmm 

TtVU 

Hiblll 

wfssm 

.086 

FAiUfTmO—  1 

wmm 

11W3 

■fcWlU 

mn&n 

Mtw»a 

massm 

Overall  Average 

.142* 

.148* 

.152* 

.160* 

.166* 

.179* 

.180* 

.161* 

AV.  Perpendicular 

.071 

.174 

.078 

.080 

.083 

.087 

.090 

.080 

^values  estimated  from  the  value  for  Perpendicular 
orientation 


TABLE  VIII 


Temperature:  -16.5°C 

Q 

Frequency  (xlO  Hz) 


28 

30 

32 

34 

36 

38 

RSI 

Average 

1  1 

.054 

.054 

.054 

.053 

.053 

W»m 

.053 

.053 

O 

2.85 

EH 

.035 

.037 

.038 

.039 

.041 

.042 

.044 

.039 

-L 

imcn 

PBEB1 

BSH1 

HQVE.Hh 

WEBM 

HilcMi 

ianai 

RH 

U 

i  1 

.053 

.054 

.054 

.055  1 

.056 

■rota 

EH 

•H 

c 

3.40 

EH 

.039 

.04J 

.040 

.041 

.041 

.042 

.042 

H 

•H 

rH 

m 

\wmm 

msam 

.027 

wcssm 

HW.1 

■ftw m 

.027 

in 

1 1 

.130 

.110 

.104 

.098 

.090 

Biumi 

7.20 

AV. 

.084 

.079 

.075 

.070 

.066 

.060 

.056 

.070 

L 

.036 

.039 

.032 

.030 

.028 

■BFU 

.030 

AV.  Parallel 

.079 

.077 

.075 

.073 

.071 

.069 

■ 

1  H 

Overall  Average 

IHiJJJi 

.052 

.051 

mm 

.050 

AV.  Perpendicular 

.027 

.027 

.027 

.028 

.028 

.028 

.028 

.028 
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Comparison  of  the  Average  Values 
of  Loss  Tangent 

TABLE  IX 


Temperature:  -21.5°C 


Frequency  (xlO^Hz) 


28  30  32  34  36  38  40  Average 


Salinity  %o 

2.85 

1  1 

.034 

.036 

.038 

.039 

.041 

.043 

.044 

.039 

AV. 

.023 

.025 

.027 

.029 

.031 

.033 

.035 

.029 

-L 

\wmmm 

.013 

■mm 

.018 

■I&ll 

wmm 

tfjFPM 

•  018 

3.40 

n  " 

wmm 

mum 

.036 

.036 

■RiliS 

wmm 

WffEl 

.027 

wmm 

.026 

± 

.022 

.021 

HIVJII 

Hoil 

BB1 

.016 

.019 

7.20 

1 1 

.070 

.068 

.067 

.065 

.064 

.062 

.060 

.065 

AV. 

.048 

.046 

.045 

.043 

.034 

.043 

L 

\mwmm 

mmM 

WV$¥m 

R^l 

wmm 

KiiLi 

HH 

AV.  Parallel 

.047 

.047 

.047 

.047 

.047 

.047 

■TOMi 

mm 

Overall  Average 

.033 

.033 

.033 

•  C33 

.033 

.033 

.033 

.033 

AV.  Perpendicular 

.019 

.019 

.019 

.019 

.019 

wmm 

msm 

7019" 

TABLE  X 

Temperature:  -25°C 

Frequency  (xlO^Hz) 


28 

30 

32 

34 

36 

38 

40 

Average 

1  1 

.023 

■n>» 

wmm 

E221 

E229I 

wmm 

.028  [ 

2.85 

AV. 

.016 

.018 

.020 

.022 

.025 

.027 

.029 

wmam 

nHii 

.011 

wmm 

BIU 

■Void 

.022 

■EM 

7017 

>* 

1 1  - 

.027 

wmm 

.027 

.028 

.028 

.028 

.028 

.028 

•H 

3.40 

AV. 

.024 

.024 

.024 

.023 

.023 

.023 

.023 

.024 

•H 

-L 

KRil 

mrw 

.018 

.018 

.019 

i— 1 

CO 

1  1 

.045 

.044 

.052 

.055 

.059 

.062 

.065 

.055 

7.20 

RHI 

.033 

.034 

wmm 

msm 

.039 

.041 

.033 

L 

.020 

.019 

.019 

.018 

.017 

.017 

.016 

.018 

AV.  Parallel 

.032 

wmm 

■TOM 

K221 

mm 

wmm 

.037 

Overall  Average 

.024 

.025 

.026 

.028 

.029 

.030 

.031 

.028 

AV.  Perpendicualr 

.016 

.017 

.017 

.018 

.018 

.019 

.019 

B!2£aH 
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Comparison  of  the  Average  Values 
of  Loss  Tangent 

TABLE  XI 


Temperature:  -32°C 

Frequency  (xlO^Hz) 


28 

30 

32 

34 

36 

38 

40 

Average 

1 1 

.002 

.003 

.004 

.006 

.006 

.007 

BiliM 

.005 

2.85 

AV. 

IHHH 

■ 

■Hi 

mm\ 

o 

L 

irniiyq 

.002* 

.002* 

.003* 

.003* 

.004* 

004* 

■ii  m\ 

.OUT- 

TiTiM 

.009 

.009 

.010 

.010 

.009 

4J 

•H 

3.40 

AV. 

.006 

.007 

.007 

.007 

.007 

.007 

.007 

.007 

a 

♦H 

L 

.006 

mmm 

.006 

hiim 

.005 

RS9 

Hun 

tmm\ 

r*H 

■ll  ■[ 

.009 

.009 

.010 

.010 

.011 

.011 

.012 

.009 

(/) 

7.20 

AV.  | 

.006 

.006 

.006 

.007 

.007 

.008 

.008 

.007 

_ 

pm 

.003 

.003 

.004 

.004 

HlIlM 

H’l’La 

fCTiJH 

AV. 

\mmm 

wmm 

mm 

.008 

.009 

.009 

.010 

W»TiT:— 

Overall  Average 

itiPM 

.005 

.006 

.006 

.007 

.007 

.006 

AV,  Perpendicular 

iii'f.M 

.004 

.004 

.004 

.004 

.004 

.004 

.004 

*estimated  from  values  for 

Parali 

.el  Or: 

.ent 

TABLE  XII 

Tan  6  vs  Frequency 


f 

— 

_L 

■ 

_L 

■ 

_L 

28 

.058 

.076 

.020 

.081 

.029 

.111 

.032 

30 

.060 

.080 

.023 

.003 

.029 

.113 

.032 

32 

.062 

.083 

.026 

.087 

.030 

.116 

.032 

34 

.064 

.085 

.028 

.090 

.029 

.118 

.032 

36 

.065 

.088 

.030 

.092 

.030 

.120 

.032 

38 

.067 

a 

.032 

.095 

.031 

.121 

.032 

40 

.069 

.034 

.098 

.031 

.123 

M 
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3.70 


Figure  11.  Change  in  the  Real  Part  of  the  Complex 

Permittivity  With  Temperature  and  Frequency 
Averaged  for  all  Salinities. 
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Figure  12.  Change  in  the  Real  Part  of  the  Complex 

Permittivity  with  Temperature  and  Salinity 
for  a  Constant  Frequency  of  34  Ghz. 
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Figure  13.  C.ange  In  the  Real  Part  of  the  Complex 
Permittivity  With  Frequency  and  Salinity 
Averaged  for  all  Temperatures. 
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Frequency  (  x  10  Hz) 


Figure  14.  Change  in  Loss  Tangent  with  Frequency,  Salinity,  and 
Orientation} Averaged  for  all  Temperatures. 
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Figure  15.  Change  In  Loss  Tangent  with  Salinity, Frequency »and 
Orientation;  Averaged  for  all  Temperatures. 
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Figure  16.  Change  in  Loss  Tangent  with  Temperature  and  Orientation 
Averaged  for  Frequency  and  Salinity. 
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Figure  17. 

Change  in  Loss  Tangent  with 
Temperature, Orientation,  and 
Frequency;  Averaged  for  all 
Salinities. 


VI.  DISCUSS  ION 


An  examination  of  the  grapnic  presentation  of  the  data  for  the 
real  part  of  the  dielectric  constant  (  e')  shows  that,  though  differ¬ 
ences  do  exist  for  individual  conditions,  the  effect  of  brine  inclu¬ 
sion  orientation  is  small  and  not  detectable  within  the  accuracy  of 
these  experiments.  It  has  been  stated  that  sea  ice  can  be  considered 
as  a  mixture  of  two  dielectrics,  ice  and  brine,  and  that  the  complex 
permittivity  will  be  a  function  of  the  interaction  of  these  two  com¬ 
ponents.  At  the  warmest  temperature,  -7°C,  and  highest  salinity, 

7.2  °/oo,  covered  by  these  experiments  the  brine  value  was  approximately 
6%  of  the  total  volume  of  the  sea  ice  (see  Figure  2) ,  compared  to 
2.5%  for  2.85  °/oo  sea  ice  at  -7°C.  The  percent  brine  by  volume  drops 
to  one  tenth  of  these  values  at  -32°C.  2  Figure  11  shows  that  the 
average  value  to  the  real  part  of  the  complex  permittivity  decreases 
with  decrease  in  temperature,  for  the  higher  frequencies  measured. 

Figure  12  shows  that  this  effect  is  due  primarily  to  the  samples 
measured  at  7.2  °/oo  salinity.  It  should  be  noted  that  the  errors 
in  these  measurements  are  independent  of  salinity  and  frequency. 

Further,  the  values  represented  in  Figure  12  for  34  GHz  also  re¬ 
present  the  averages  for  all  frequencies.  Figure  13  shows  the 
overall  increase  of  the  real  part  with  increase  in  frequency.  These 
effects  will  have  to  be  examined  in  more  detail  before  conclusions 
can  be  drawn.  No  physical  explanation  for  this  observation  has  been 
found. 
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The  most  interesting  result  from  this  series  of  measurements  is 
seen  in  the  marked  effect  of  brine  inclusion  orientation  on  the  loss 
tangent  in  sea  ice.  An  examination  of  the  data  presented  in  Figures 
A15  through  A30  of  the  Appendix  shows  that  without  exception,  the  par¬ 
allel  orientation  of  brine  inclusions  with  respect  to  the  electric 
field  in  the  wave  guide  caused  the  loss  tangent  to  be  greater  than 
perpendicular  orientation.  It  has  been  shown  by  a  number  of  investi¬ 
gators  that  at  frequencies  above  approximately  1  x  10Q  Hz  dielectric 
losses  due  to  ionic  conduction  are  negligible  and  loss  can  be  princi¬ 
pally  associated  with  dipole  rotation  of  the  molecules  of  the  dielectric 
material.12,  1  3<  114  In  a  crystalline  substance  such  as  pure  ice,  mole¬ 
cular  rotation  is  restricted  by  the  nature  of  the  molecular  bonding 
and  dielectric  loss  is  small.  In  liquid  water,  dipole  rotation  is 
significant  and  causes  dielectric  loss  to  be  high  in  the  frequencies 
of  interest  in  this  experiment.  1  2-  1 3i  14,  It  has  been  shown  that  the 
shape  of  the  included  medium  in  mixtures  of  dielectrics  has  an  effect 
on  the  dielectric  properties  of  the  mixture.14, 15,  This  theory  should 
be  examined  further  for  the  specific  case  of  sea  ice. 

Figure  14  shows  that  the  loss  tangent  of  sea  ice  increases  with 
frequency  and  that  there  is  a  pronounced  difference  between  the  differ¬ 
ent  orientations  for  all  salinities  measured.  Figure  15  shows  that 
loss  increased  with  salinity  and  again  shows  the  increase  with  fre¬ 
quency.  Figure  16  shows  the  large  change  in  loss  with  temperature 
and  orientation.  Figure  17  shows  (his  change  with  frequency  as  an 
added  parameter. 


Figure  18  shows  a  plot  of  ln(  tan  5)  as  a  function  of  1/T, 
temperature  being  measured  in  degrees  Kelvin.  It  would  seem  appar¬ 
ent  from  this  presentation  that  a  dramatic  change  occurs  within  the 
sea  ice  in  the  temperature  range  from  -21.5°C  to  -25°C,  since  the 
slope  of  the  straight  line  indicates  the  activation  energy  of  the 
physical  process  which  results  in  loss. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


We  can  conclude  that  the  average  value  of  the  real  part  of  the 
relative  complex  permittivity  for  this  band  of  frequencies  is  3.08 
-  0.06.  This  compares  favorably  with  limiting  values  for  pure  ice 
at  high  frequencies  calculated  by  Auty  and  Cole  to  b'?  3.10  and  by 
Lamb  to  be  3.17. 12  These  values  were  extrapolated  from  experiments 
done  at  lower  frequencies.  Pure  ice  should  be  examined  at  the 
frequencies  of  this  experiment  for  comparison.  The  increase  in  the 
real  part  seen  in  these  measurements  with  the  higher  frequencies 
at  low  temperatures  should  be  examined  more  closely. 

These  measurements  show  strong  indications  that  the  loss  tangent 
of  sea  ice  is  a  function  of  the  orientation  of  the  brine  cell  vertical 
axis.  The  physical  process  behind  this  effect  should  be  examined  more 
closely. 

Figure  18  represents  a  plot  of  the  expression 

In  (tan  5 )  =  InA  -  ^  (24) 

KT 

where  -  -j^-  represents  the  slope  of  the  curve  and  can  be  related  to 
the  activation  energy  of  the  loss  process  in  the  sea  ice.  Calculations 
show  that  the  average  slope  for  all  salinities  between  -7°C  and  -25°C 
is  -5800.  However,  the  average  slope  between  -25°C  and  -32°C  is  -13,400. 
There  is  an  indication  that  a  change  occurs  in  sea  ice  between  -21°C 
and  -25°C  which  affects  the  nature  of  the  loss  process.  This  effect 
should  be  investigated  further.  Empirical  expressions  for  the  loss 
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tangent  of  sea  ice  at  approximately  5.0  °/oo  can  be  derived  from  the 
information  in  Figure  18,  for  the  two  temperature  ranges.  For  temper¬ 
atures  above  -259C  this  is 

In  (tan  6)  *  20.13  -  5800  |  (25) 

For  temperatures  below  -25SC  this  expression  will  be 

In  (tan  6)  -  50.43  -  13400  -  (26) 

T 

Temperature  in  these  equations  is  expressed  in  degrees  Kelvin.  These 
expressions  will  need  refinement  for  increased  accuracy,  but  they 
should  provide  a  good  first  approximation  of  the  loss  factor  for  sea 
ice  below  approximately  -5°C  and  in  the  vicinity  of  5  °/oo. 
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APPENDIX 


TABLE  AI  Data  for  Sea  Ice  of  Salinity  2.85% 

TABLE  All  Data  for  Sea  Ice  of  Salinity  3.40% 

TABLE  AIII  Data  for  Sea  Ice  of  Salinity  7.20% 

FIGURES  Al 

to  A15  Least  Squares,  Best  Straight  Line  Fit  for  the  Real 

Part  of  the  Complex  Permittivity,  All  Samples 

FIGURES  A16 

to  A30  Least  Squares,  Best  Straight  Line  Fit  for  the  Loss 
Tangent,  All  Samples 


TABLE  AI.  DATA  POINTS  FOR  SEA  ICE 
OF  SALINITY  2.85  °/oo 


Temperature  Frequency  (x  109  Hz) 


perpendicular 

parallel 

perpendicular 

parallel 

-7#C 

30.02 

3.07 

.074 

33.58 

3.47 

.103 

33.56 

3.13 

.049 

37.18 

3.09 

.083 

29.57 

3.22 

- 

29.60 

2.20 

- 

32.21 

3.43 

- 

32.30 

2.88 

- 

34.54 

4.13 

- 

35.06 

3.29 

- 

37.90 

2.87 

- 

-16.5#C 

29.23 

2.81 

.017 

30.62 

2.92 

.026 

34.15 

2.93 

.024 

34.28 

3.19 

.032 

37.98 

3.01 

.033 

28.99 

3.27 

- 

30.16 

2.83 

.048 

30.92 

3.07 

.048 

31.34 

3.02 

.069 

32.68 

3.26 

.070 

33.81 

2.77 

- 

35.23 

2.91 

.063 

36.24 

3.42 

.069 

37.93 

2.58 

.052 

38.20 

2.79 

.044 

38.67 

2.89 

- 

-21. 5°C 

29.28 

2.76 

1 

30.75 

3.00 

1 

34.25 

2.82 

.013 

34.43 

3.06 

38.11 

3.22 

.029 

29.22 

2.36 

.037 

30.00 

3.05 

.037 

30.20 

2.81 

.034 

31.58 

4.13 

.044 

32.77 

3.05 

.041 

32.80 

3.00 

.036 

(continued  on  the  next  page) 


Table  AI  (continued):  Data  Points  for  Sea  Ice  of  Salinity  2.85  °/oo 


Temperature  Frequency  (xlO^  Hz) 
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TABLE  All:  Data  Points  for  Sea 

Ice  of  Salinity  3.40  °/oo 


Temperature  Frequency  CxlO9  Hz 


Tan  6 


p«rp.  I  Por 


TABLE  All:  Data  Points  for  Sea  Ice  of  Salinity  3. AO  °/oo 
(continued) 
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TABLE  All I: 


Data  Points  for  Sea  Ice  of  Salinity  7.20  °/oo 


(continued) 
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. - Perpendicular  Orientation 

' - Parallel  Orientation 

- Average 


- Perpendicular  Orientation 

Salinity:  2.85%o  Key: - Parallel  Orientation 

Temperature:  -16.5°C  - Average 


Figure  A-2: 


z'  versus  Frequency 


58' 


i 


Salinity:  3.4X0  - Perpendicular  Orientation 

Temperature:  -16.5°C  Key: - Parallel  Orientation 

-  Average 


Figure  A- 10:  e'  vs.  Frequency 


Salinity:  7.2% 
Temperature:  -7°C 


Figure  A-ll:  c'  versus  Frequency 


-67- 


•v 

t  v 


Figure  A-12:  e'  versus  Frequency 


68' 


Frequency  (x  10  Hz) 


Salinity:  7 .2%a 

Temperature:  -21.5°C 


Unity 


Frequency  (  x  10  Hz) 


Salinity:  2.S5% 
Temperature:  -7°C 


Frequency  (x  10  Hz) 


Frequency  (  x  10s  Hz) 
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Figure  A-18:  Tan  6  versus  Frequency 
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Salinity:  3.40% 
Temperature  -21.5°C 
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Temperature  :  -16.5°C 
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